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Abstract 

The exact solutions for linear cosmological perturbations which have been 
obtained for collisionless relativistic matter within thermal field theory are 
extended to a self-interacting case. The two-loop contributions of scalar Xcj) 
theory to the thermal graviton self-energy are evaluated, which give the O(A) 
corrections in the perturbation equations. The changes are found to be per- 
turbative on scales comparable to or larger than the Hubble horizon, but the 
determination of the large-time damping behavior of subhorizon perturbations 
requires a resummation of thermally induced masses. 
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The theory of linear perturbations of otherwise homogeneous and isotropic cosmolog- 
ical models flj plays a central role in the problem of large-scale structure formation in 
the early universe. The conventional approach to their study in cosmological models in- 
volving more complicated forms of matter than a perfect fluid [p]] is based on the coupled 
Einstein-Boltzmann equations [0-|J, which are usually solved numerically In Ref. |7J, a 
novel framework has been developed which employs perturbative thermal field theory || 
to derive self-consistent and automatically gauge- invariant perturbation equations. In this 
formalism, the connection of the perturbed metric and the perturbed energy-momentum 
tensor is provided by the thermal graviton self-energy. At one-loop order, the leading high- 
temperature contributions to the latter have been obtained in Ref. (see also Ref. |TIJ) in 



momentum space, and conformal invariance allows one to transform these results directly 
to a curved space with vanishing Weyl tensor. The one-loop high-temperature contributions 
describe relativistic collisionless matter interacting with the metric perturbations, and the 
resulting perturbation equations have been shown [] f to be equivalent to a gauge-invariant 



reformulation of the Einstein- Vlasov equations ||12|| . Moreover, analytic solutions for scalar, 
vector, and tensor perturbations have been found for a spatially flat Friedmann-Robertson- 
Walker (or, Einstein-de Sitter) cosmological model in the general case of a two-component 
system containing also a perfect fluid part [f?|, [TT]JT5| ] , as well as for higher-dimensional models 

0- 



In this Letter, we investigate the effects of weak self-interactions on the exact solutions 
that were found for relativistic collisionless matter. For simplicity we consider scalar particles 
with A0 4 interaction Lagrangian. The leading self-interaction effects are described by two- 
loop corrections to the thermal graviton self-energy in the high-temperature limit. The 
corrections are such that the ensuing perturbation equations can still be solved exactly in 
terms of rapidly converging power series. 

The background space-time of the Einstein-de Sitter cosmological model, which we shall 
consider, is given through the line element 

ds 2 = S 2 (r) (-dr 2 + 5ijdx l dx J ) , (1) 

with r being the conformal time which measures the size of the horizon in comoving co- 
ordinates, Rh = t. A thermal distribution of relativistic (i.e., effectively massless) scalar 
particles with A0 4 self-interactions gives rise to a background energy- momentum tensor 

f\ = u»u u {E + P) + F<f£ , u, = S6° , (2) 

with E = 3P. Including order A, T^ v is determined by the diagrams of Fig. 1, which at 
finite temperature T yield 



i{s - '> = (s - b) t> + 0(AS/2) • 



(3) 



The scale dependence of the mean energy density is given by E(S) = E(S = l)^ -4 , with 

S OC T. 

Self-consistent perturbations of the metric = + 5g^ u are found by varying the 
Einstein equations, 
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with 



bGT = 5(R txv - -g^R) = SttGST^ (4) 



5T ^ X) = I d4y d^ 6ga " iy) - (5) 



Hence, ST^ U is determined by the graviton self-energy 



where V contains all contributions to the effective action besides the classical Einstein-Hilbert 
action. 

T is a nonlocal functional of the metric. With a flat background metric, the graviton 
self-energy H^ap is most easily evaluated through momentum-space Feynman rules. In 
the high-temperature limit, V is conformally invariant |7] which makes it possible to directly 
transform the flat-space results to conformally flat cosmological models such as the Einstein- 
de Sitter one of Eq. ([!]) according to 

= s ~ 2(r) / ^y etk{x ~ x ' ] n^^L ^ 2 ( r ') • ( ? ) 

The one-loop contribution to f[^ upcr has been evaluated in Ref. at temperature 
T ^> k , k, which is relevant for the case of cosmological perturbations [the momentum 
scale is set by the inverse Hubble radius ~ (r/mpi anc k)T]. It describes collisionless mat- 
ter, interacting solely with the metric perturbations. Since collisionless matter can sustain 
anisotropic pressure, this medium admits a rich spectrum of possible cosmological solutions. 
Those have been studied in great detail in Ref. [|ll|] by two of the present authors. 

In order to investigate the effects of weak self-interactions, one has to include higher loop 
corrections. With A0 4 interactions, the order A contributions are contained in the two-loop 
diagrams shown in Fig. 2. Evaluating again the high-temperature limit T ^> ho, k (full 
details will be presented elsewhere) gives a contribution to n Miypcr (fc) that again satisfies the 
Ward identity corresponding to conformal invariance of T which is necessary for the simple 
transformation law of Eq. (|7|). 

Since we did not have to single out a particular gauge, the perturbed Einstein equations 
(|) involve only gauge invariant combinations of the perturbed metric components and can 
be given entirely in terms of Bardeen's gauge invariant potentials JT5 . Expanding the latter 



in Fourier modes (plane waves with respect to the conformally flat background of Eq. (TI])) 
then leads to equations depending only on the conformal time r, which can be combined 
with the wavelength of the metric perturbations £ = 2%/k into the dimensionless variable 

*^r=?fS. (8) 

x/n = 1 defines the point in conformal time when half a wavelength fits inside the growing 
Hubble horizon. 
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Because of the nonlocality of the effective action T, the resulting equations are integro- 
differential equations rather than ordinary differential equations, involving a convolution 
integral of n M ^ p<J (x, x') with the metric perturbations. 

In the case of scalar metric perturbations, there are two independent equations for the 
two gauge invariant potentials and $a of Bardeen [I5j . Using instead the combinations 

1 



H, 



n 



(9) 



which are related to the gauge-invariant (cf. |I5| ) density contrast e m = x 2 $/3 and to the 
anisotropic pressure 7r anis = x 2 U, one equation (the 00-component of Eq. (B)) is of the form 



x 2 - 3)<$>{x) + 3x$'(x) - 6Tl(x) = 
-12 J dx'K(x - a/)($ + Tl)'(x') 



(10) 



while the second is an ordinary differential equation by virtue of the vanishing of the trace 
of the energy-momentum tensor in our ultrarelativistic setting, 



4 12 2 

$" + + _$ + _n' - -n 

x 3 x 3 



. 



All the remaining components of the perturbed Einstein equations 
binations of Eqs. (|l(J) and flTTl) . 



(11) 

turn out to be com- 



The one-loop contribution to the kernel K is determined by the discontinuity of the 
one-loop graviton self-energy Tl 0000 (k , k) in Eq. (0). It reads @ (u = ko/k) 



K x (x 



x 



i 

2^ 



to 



-oo+ie 



j (x - x')6(x - x'), 



(12) 



with jo(x) = sin(x)/a;. The step function gives an upper bound of x in the integral in 
Eq. (PUD; a lower bound at some initial x can be adopted by adding an inhomogeneous 
term Y^=olnK^ n \x — Xq) to the convolution integral in Eq. ([T0|), encoding the infinite 
number of initial conditions one has to impose by specifying the likewise infinite number of 
moments of the initial particle distribution JTTJ . 

The two-loop contributions have more complicated structure, but can still be evaluated 
exactly, 



Ko(x — x' 



5i\ 
32tt 3 



oo+ie 



dcu e"^-*') 



x 



oo \n 2 



OJ 



In 



u + 1 



2oo 



CO 



CO 



CO* 



5A 

8^ 



{2k — j — cos)(x — x')8(x — x) 



(13) 



where we introduced 
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k(x) 



- [sin(x)Si(2x) + cos(x){Ci(2x) - 7 - ln(2x)}] 



2E 



m=0 



(2m + 2)! 



■E 

3=0 



2j + l 



(14) 



with Si and Ci being the sine and cosine integral [16]. In Eq. (|I~3"D, there are two contributions 
which differ in form from the one-loop kernel: the first term in the integrand, which arises 
exclusively from the last diagram in Fig. 2, and the last term, which is due solely to the 
third diagram. 

With the explicit result for K 2 , which has a well-behaved power series representation, 
the same methods as in Ref . [0,0] can be used to solve the coupled equations (|T0f) and ( |i~l~l) . 
Since K 2 (x) oc x 4 + 0(x 6 ), K = K\ + K 2 still satisfies 



K{0, 



K(0+) + 3K"(0+) = 



(15) 



as did K\. Thereby the conditions [JTTJ] on the initial data (summarized by the set of 7 n ) 
following from demanding regularity at x = xq = are left unchanged. 

In addition to regular solutions, there exist also ones that are singular as x — > 0. In 
the collisionless case those behave as $, II ~ x~ 5 / 2 cos[ 
tions even on superhorizon scales [11,17,18 
interactions, to wit, 



'27/20 ln(x)], thus showing oscilla- 
This behavior is modified by the weak self- 



~ x 5 / 2 cosf 



I-^ ln(x) 



7T- 



(16) 



up to a constant shift of phase. 

In Ref. |7J , it was pointed out that the coupled system of Eqs. ( [TOD and ( [TT]) can be solved 
exactly for xq = 0, yielding power series representations with infinite radius of convergence. 
This still holds true after adding the self- interaction term K 2 . For the particular initial 
condition determined by 7„ = 0, n > 1, the first few terms in the regular solutions are given 
by 



$(a;) 



/28- IOA/tt 2 
70 [ 5^ 



(17) 



4644 - 3735A/vr 2 + 350(A/ 



7T 



2\2 



315(54 - 5A/tt 2 ) 



X 



± 



U(x) 



To 



4-5A/ 



7T 



2808 - 5445A/7T 2 + 700(A/tt 2 ) 2 2 



630(54 - 5A/tt 2 



x =F 



With A = 1, which corresponds to a prefactor rs 0.06 in the last line of Eq. ([T3|), the 
density contrast e m = x 2 $/3 is plotted in Fig. 3 by the dashed line and compared with the 
collisionless case (full line). On superhorizon scales, 1 C 1, the solutions are dominated by 
the constant modes in $ and II. The amount of anisotropic pressure associated with a given 
density contrast is determined by II/$, which is seen to be reduced with increasing A, in 
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agreement with intuitive expectation, since the collision-dominated case of a perfect fluid 
has II = 0. 

For large i>1, the asymptotic behavior of the density contrast e m (x) is proportional 
to the one of K(x). In the collisionless case K(x) = K\{x) ~ sin(x)/x, whereas in the other 
extreme case of a perfect fluid e m ~ sin(x/v^3). With weak self-interactions, the asymptotic 
behavior is indeed modified towards weaker damping (see Fig. 3 and 4), which ceases com- 
pletely for sufficiently large x where K(x) ~ Acos(x). However, there the correction K-i is 
overtaking the lowest order term Ki, which signals a possible breakdown of perturbation 
theory. 

Inspecting the analytic structure that gives rise to the different asymptotic behavior of 
K2, one finds that the discontinuity in the integrand of Eq. (|i~3"D is singular on the light- 
cone u = ±1. There is a logarithmically singular term which is responsible for k'(x) ~ 
sin(x) ha.(x)/x, but there are even poles at u — ±1, which contribute the term involving 
the cos(a;). However, the scalar particles, whose rest masses are negligible at sufficiently 
high temperature, acquire thermal masses m = \f\T . A resummation of these those will 
remove the singularities at u = ±1, which are associated with the masslessness of the scalar 
particles. 

The logarithmically singular integrand of Eq. ( |T3"D will be regularised at uj = ±1, but 
since the logarithmic singularity is integrable anyway, the resummed result will not change 
dramatically; the poles at uj = ±1, however, disappear completely after a resummation of 
the thermal masses. 

Adding a thermal mass term to the Lagrangian and subtracting it again as a higher-order 
counter-term has the effect of replacing all the propagators in Fig. 2 by massive ones and 
subtracting out the third and fourth diagram, which were precisely the ones contributing 
the pole term to Eq. ([T5]) . This resummation of mass insertions replaces the leading-order 
result Ki(x) — jo(x)9(x) by 

(18) 




up to terms whose amplitude is down by factors of A. With a non-zero m = y/XT, the 
integrand is now seen to vanish at to = ±1, whereas before it was a constant; in fact, it 
vanishes together with all its derivatives, but rapidly recovers the bare one-loop value away 
from the light-cone. This is a negligible effect for small x, but the large x behavior is 
dominated by the (non-analytic) integration region |u;| ~ 1. 

Eq. ( |T%| ) can be evaluated by a Mellin transform [T^] which yields 



rw = «w|^x (19) 



f Afc/2 (-l) 1+fc (fe-4)C(fc-3) (2^ 
where in the term with k = 4 one has to substitute (k — 4:)((k — 3) —>■ 1. For small x, 
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Krix) = 6{x) (j Q (x) - cos(x) +0(A 3 / 2 ) (20) 



is a good approximation; for i > 1, on the other hand, the complete function K[ es '(x) turns 
out to decay even slightly faster than jo(x), oscillating with a reduced phase velocity 

' :>A 0(A 3 / 2 ). (21) 



In 2 



A first approximation to a fully resummed calculation is therefore to change 



5A 

Ki(x) oc j {x) -> j ((l - ^) x ) 



which reduces 



K 2 {x) -> —(k - j ){x)6{x) 



5X 
4^ 

up to higher orders in A. However, the latter of the relations in Eq. fll5|) is now violated 
at the higher order 0(A 2 ), which drastically restricts the solvability of the perturbation 
equations, so further modifications are needed. A natural solution is to change the phase 
velocity in K 2 by the same amount as in K\ and to adjust the prefactor in K 2 appropriately. 
This leads in a unique manner to the remarkably simple result 

K™-{x) = (j + - Jo)) (ike) 9{x). (22) 

Up to and including order A, this is the same as K\ + K 2 . But now the asymptotic behavior 
of e m (x) for large x, which is linked to K(x), is ~ k'(vx) ~ sin(fx) ln(x)/x. Thus there is 
still strong damping of subhorizon-scale oscillations, which is only logarithmically weaker 
than in the collisionless case. Moreover, the phase velocity v is now smaller than 1, as one 
may expect from the comparison with the strongly collisional case of a perfect radiation 
fluid, where (undamped) acoustic waves oscillate with v = l/y/3. In Fig. 4 the effects of the 
above resummation are shown for x/tt > 1. 

Similar issues arise also in the case of vector (rotational) as well as tensor perturbations 
(primordial gravitational waves). For small x, weak self-interactions affect mainly the sin- 
gular solutions, and the same shift of frequency of the superhorizon oscillation occurs as 
in Eq. (|16D ; for large x, vector and tensor perturbations differ from the scalar case in that 
their asymptotic behavior is not changed as much. Subhorizon-scale vector perturbations 
still have vorticity decaying ~ sin(x)/x, and resummation changes only the phase velocity; 
tensor perturbations are stable against the kind of resummation performed above, and have 
the same asymptotic behavior as the collisionless ones, which indeed is already the same as 
in the perfect fluid case. All these findings as well as a generalisation to more-component 
systems will be presented in detail in a forthcoming publication. 
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FIGURES 




FIG. 1. One- and two-loop diagrams for the energy-momentum tensor. Wavy lines denote 
external gravitons and straight lines scalar particles. 




FIG. 2. The one- and two-loop diagrams for the graviton self-energy up to and including 0(A). 
Only the last diagram contributes to the At' -term in K%. 
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x/pi 



FIG. 4. In addition to the lines in Fig. 3, the dotted line shows the solution with the resummed 
kernel (A = 1). 
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